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Abstract 

We study the behavior of Heisenberg, antiferromagnetic, integer-spin chains in the pres- 
ence of a magnetic field exceeding the attendant spin gap. For temperatures much smaller 
than the gap, the spin chains exhibit Luttinger liquid behavior. We compute exactly both 
the corresponding Luttinger parameter and the Fermi velocity as a function of magnetic 
field. This enables the computation of a number of correlators from which we derive the 
spin conductance, the expected form of the dynamic structure factor relevant to inelastic 
neutron scattering experiments, and NMR relaxation rates. We also comment upon the 
robustness of the magnetically induced gapless phase both to finite temperature and finite 
couplings between neighbouring chains. 

1 Introduction 

The existence of a gap in one-dimensional, integer-spin, Heisenberg antiferromagnets was first 
predicted by Haldane [T]. He found that such spin chains can be mapped onto a gapped field 
theory in the large-spin, continuum limit. A variety of checks imply that this behavior persists 
down to spin s = 1. A spin-1 chain with a specific {S ■ S)"^ coupling has been rigorously shown to 
exhibit a spin gap |4 . While at a differing value of the {S ■ S)"^ coupling, the spin chain is gapless 
IS] , this critical point is believed to be unstable in the two-parameter space of couplings. Gapless 
behavior thus only arises as a product of fine-tuning. Numerous numerical studies carried out on 
spin-1 chains observe a gap [HI. Experimentally, inelastic neutron scattering studies on a number 
of quasi-one-dimensional spin-l-chain materials are consistent with a finite spin gap [ZlElin]- 

The physics underlying the gap is particularly robust: related systems such as two-leg spin- 
1/2 or Hubbard ladders also exhibit a gap to spin excitations ilOj. Roughly speaking, integer-spin 
composites form across the rungs of the ladder making it into an effective integer-spin chain. Both 
the ability to fabricate these materials and their relationship to high cuprate superconductors 
have made them the focus of intense theoretical and experimental study jTH I12j. 
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The centrality of the gap is no less when these materials are subjected to an applied Zeeman 
field, H . In a field exceeding some critical value, Hf. ~ A, the excitation spectrum of integer- 
spin chains changes dramatically and the gap vanishes ^3^^. At i7 = the ground state is a 
singlet. The elementary excitations above the ground state are three massive spin-one bosonic 
magnon modes. In contrast, when H exceeds He-, the gap of one of the magnons closes and the 
ground state of the spin chain begins to fill in with gapless excitations. If interactions between 
the gapless magnons were completely absent, the excitations would collapse into a condensate 
of free bosons. On the other hand if interactions were perfectly repulsive at any energy scale, 
the bosons could be thought to possess a hard core and so form a gas of free fermions with some 
Fermi surface. But in fact the interactions are expected to lie midway between these extremes 
and what one ends up with is equivalent to an interacting sea of fermions. This of course is 
nothing more than a Luttinger liquid. The purpose of this article is to study its properties. 

The field theory describing the continuum limit of integer-spin chains is the 0(3) non-linear 
sigma model (NLSM) without topological term ^ |2]. The 0(3) NLSM also describes half- 
integer-spin chains but in this case the topological term is present with coupling 6 = n. The 
presence of the topological term leads to gapless behavior. Indeed at low energies the 0(3) 
NLSM with 9 = -K fiows to the critical theory, SU {2)i (for a review see |2I). 

Unlike its ancestral theory, the Heisenberg spin-1 chain, the 0(3) sigma model has the virtue 
of being integrable |Sl E] • While aspects of the physics of integer-spin chains have been stud- 
ied using the integrability of the 0(3) NLSM[l31 El El E|; spin chains remain incompletely 
understood. While thermodynamic properties of an integrable model are generically accessible, 
correlation functions are not. There do exist a variety of techniques to compute correlators. At 
zero temperature, truncated 'form-factor' computations are able to access exact information of 
the low-energy properties of the spectral functions ^3 E] ■ These techniques have recently been 
extended to compute exact low temperature expansions of these same quantities jI9j. There also 
exists an elegant semi-classical approach 1^ predicated upon combining ultra-low-energy 
information from the quantum theory with older approaches of computing correlators in classi- 
cal systems [221 • However, all of these techniques require the energy scales in the problem (in 
particular, the temperature, T, and the applied field, H) to be much smaller than the gap, A. 

Fortuitously there exists an alternate approach which allows the computation of correlation 
functions when H > He- The control over the thermodynamic properties of the model that 
integrability affords allows the computation of the Luttinger parameter. While the Luttinger 
parameter characterizes the strength of the interactions between excitations, it is a far more 
central quantity in that it provides a near complete description of the low-energy structure of 
the theory. Together with the knowledge of the Fermi velocity, a quantity also easily accessible 
with integrability, a host of information can be determined, including the computation of specific 
correlators and their scaling exponents. Such techniques have been used to study the sine-Gordon 
model in the presence of a chemical potential exceeding the gap (i.e. the mass of the sine-Gordon 
solitons) 12311^ • Similarly the low-energy structure of doped generic Hubbard ladders/armchair 
carbon nanotubes has been studied in [12j utilizing this formalism. The particular derivation of 
the Luttinger parameter in this paper is based upon the treatment in .25j . 

The ability to compute the Luttinger parameter together with the Fermi velocity is predicated 
upon some generic properties of integrable models. Most importantly, the exact eigenfunctions 
of the model's fully interacting Hamiltonian are known. With this knowledge comes a well- 
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defined notion of 'particles' or elementary excitations in tlie system. Ultimately this feature is 
a consequence of the infinite number of conservation laws possessed by the integrable model. In 
particular, particle number is conserved in any collision and multi-particle S matrix elements 
factorize into products of two-particles ones. An integrable model is a superior version of a Fermi 
liquid: a particle's lifetime is infinite regardless of distance from the Fermi surface. 

In order to appreciate these features of the 0(3) NLSM, we provide an overview of the model. 
The model is described by the action, 

S = — [ dxdtidf'nd^n), (1) 
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where n = (n^., Uy, riz) is a bosonic vector field constrained to live on the unit sphere, n ■ n = 1. 
This action is arrived at from the Hamiltonian of the spin chain, 

H = JJ2s^■S,+,. (2) 

i 

In the continuum, large-s, limit, the spin operator. Si, is related to the field, n, via 

S, = {-iysn, + M,. (3) 

n(x, t) is the sub-lattice or Neel order parameter while M describes the uniform (i.e. wavevector 
k 0) magnetization. M is related to n via 

M = -nx dtn, (4) 

and so is given in terms of the momentum conjugate to n. 

The triplet of bosons which form the low-energy excitations of the 0(3) NLSM have a rela- 
tivistic dispersion relation given by 

E(j9) = (/ + A2)V2. (5) 

Here A is the energy gap or mass of the bosons and is related to the bare coupling, g, via 
A ~ Je~^/^ . We have set the bare spin- wave velocity, Vg = 2Js (the speed of light in this 
relativistic system) to be 1. The dispersion relation of all three bosons is identical as the model 
has a global 0(3) symmetry. We stress that this relativistic invariance is a natural feature of 
the low energy structure of the spin chain. 

Our approach to the Luttinger liquid phase of the Heisenberg spin-1 chain shares some 
similarities with others taken in the literature. The NLSM has been used previously to study 
the field-induced gapless phase of spin-1 chains, albeit in the context of large- techniques j2H], 
where the 0(3) NLSM is replaced with its 0{N) counterpart. (Parenthetically, we point out 
that the large N approximation of the 0(3) NLSM have been criticized by In this work 
the authors note that the ultra low energy limit of the scattering in the 0(3) NLSM differs from 
that predicted by a large N expansion of the exact 0{N) S-matrix.) 

Beyond large-A^, the Luttinger liquid phase of the spin-1 chain has also been studied in 
the guise of analyzing spin-1/2 ladders. Bosonization techniques have been used to study the 
ladder system in a regime where the legs of the ladders are weakly coup led In the opposite 
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limit, in a regime where there are strong antiferromagnetic correlations across the ladder's rungs, 
the gapless Luttinger liquid behaviour has been explored through mapping the system onto an 
effective XXZ spin-1/2 chain In this latter case, the underlying integrability of the XXZ 
spin chain can be brought to bear upon the problem [31 . Given the simplicity of the map onto 
the spin-1/2 chain, we adapt the map to the case where there are ferromagnetic interactions 
across the rungs of the ladder (appropriate for describing a spin-1 chain). Our primary purpose 
in doing so is to come up with a more precise identification between the relevant fields in the 
theory and the bosonic degrees of freedom of the Luttinger liquid. Although we could exploit the 
integrability of the spin-1/2 chain to compute the Luttinger parameter, we do not do so. Given 
that the map to the effective spin-1/2 chain is done through a perturbative expansion where the 
importance of the higher order terms in the series is uncertain, we view extracting the Luttinger 
parameter directly from the 0(3) NLSM as more reliable. 

Haldane gap spin chains in a magnetic field have also been studied through a mapping onto 
an interacting, spinless Bose gas ^E]. This map provides a reliable description of the spin chain 
at small magnetizations. At small magnetizations, the low density of excitations forming the 
ground state interact only weakly. The scaling exponents and thermodynamic properties of the 
system are then independent of the exact nature of these interactions. Although not done by the 
authors of the integrability of the Bose gas, like the integrability of the 0(3) NLSM, could 
be used to determine the various properties of the Luttinger liquid phase. However beyond the 
low magnetization regime, the results would differ. It is interesting to note, however, that by 
fine tuning the strength of the interactions of the Bose gas, the analysis of the 0(3) NLSM, in 
particular the computation of the Luttinger parameter, can be reproduced in large degree. We 
comment upon this further in Section 3. 

In this article we take the Hamiltonian of the spin chain to be in its minimal Heisenberg form 
and so ignore (for the most part) the affects of anisotropies upon the physics. These can take 
(at least) two forms. Easy-axis anisotropies, 

AH = ^(5„)2 + Dy 5^(S,,)' + 5^(^.i)', (6) 

i i i 

of varying strengths are often found in spin-1 chain materials. Additionally, actual spin chain 
materials never take the form of an isolated chain. Rather the chains exist in three dimensional 
arrays with weak but non-zero interchain couplings, J'. Thus the chains are at best quasi-one- 
dimensional. With a finite J', there will be some correspondingly finite Neel temperature, T/v. 
Below T/v the physics will be dramatically different than described in this article. 

CsNiCls was the first material for which evidence of a Haldane gap was found 8J- This 
material suffers from the second aforementioned anisotropy with a relatively large interchain 
coupling, J' /J ~ .017. Consequently Neel order was observed to set in at T ~ 5K. A more 
promising material for the observation of a Haldane gap was found in Ni{C2H^N 2) 2^020104^ 
(NENP). For NENP, the ratio J' / J ~ 6 x 10""^ is sufficiently small that 3D Neel order has 
not been observed down to temperatures ~ 1.2ii'|ni. However this material is characterized by 
a large easy-axis anisotropy, D^/J ~ .25; -D^/A ~ 5/8 {D^ ^ Dy 0). Related materials 
Ni{C5HuN2)3iPF6) (NDMAP) and Ni{C5HuN2)2Ns{Cl0i) (NDMAZ) share similar easy-axis 
anisotropies. In terms of our analysis, these latter compounds share the additional unwanted 
feature of field-induced antiferromagnetism [32] • The Luttinger liquid that results from magnetic 
fields large enough to extinguish the Haldane gap leads to quasi-long range antiferromagnetic 
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correlations. With a small finite J', these quasi-long range correlations are promoted to full 
fledged long range order. The corresponding Neel temperature increases with applied magnetic 
fields. In a mean field framework, we compute this ordering temperature in Section 4. Thus 
at fixed temperature we expect only a finite window in the applied magnetic field in which the 
Luttinger liquid behavior will persist. 

Perhaps the material best suited to the analysis presented in this paper is AgVP2SQ. It has 
an extremely small interchain coupling, J' / J ~ 10^^ and a similarly small easy axis anisotropy, 
DJA ~ IQ-^. However it possesses a comparatively large gap, A ~ 320K. As such, high field 
measurements {H > A) have yet to be done on this material (and are unlikely to be done soon), 
opposite to the case of NENP |37j, with a much smaller average gap, A ~ 20K. 

The paper is organized as follows. In Section 2 we develop a Landau-Ginzburg description of 
the low-energy effective theory of the integer-spin chain in a magnetic field exceeding the spin 
gap. This effective theory reduces to a Luttinger liquid and so has two controlling parameters: 
the Luttinger parameter, K, and the Fermi velocity, vp- In Section 3 we show how these 
parameters can be determined as a function of the applied field, H, through consistency with 
the 0(3) NLSM. 

With this description of the low-energy theory in hand, we analyze the behavior of a number 
of correlators in Section 4. Using a Kubo formula together with our knowledge of the current- 
current correlators, we compute the spin conductance and the static susceptibility. With this 
latter quantity we have thus come full circle. We computed K and vp based upon thermodynamic 
considerations and then in turn computed correlators. From these same correlators we then 
(re) compute (consistently) thermodynamic quantities. We also study the staggered spin-spin 
correlators, quantities which would be probed both in inelastic neutron scattering experiments 
near wavevector k = n and NMR relaxation rate measurements. 

While we have already raised the issue of the effects of interchain couplings upon the stability 
of the Luttinger liquid phase, we also consider in Section 4 the robustness of this phase to 
finite temperature. Although our derivation of the effective low-energy theory is done at zero 
temperature, the conformal/scale invariant nature of the theory allows us to easily determine 
quantities at T > 0. To explore how large the temperature relative to the gap may become 
before the Luttinger picture breaks down, we study the susceptibility at finite T and H. We do 
so using a more sophisticated description of the system, the exact equations (good for arbitrary 
T and H) giving the system's free energy. 

2 Low-Energy Effective Theory for H > A 

In this section we describe the emergent Luttinger liquid behavior of a spin chain arising in 
magnetic fields larger than the gap and temperatures satisfying T <^ A. Following j^H], we 
provide a corresponding Landau-Ginzburg description. Although we work with an effective 
theory, we will be able to compute the various phenomenological parameters appearing in it by 
insisting on consistency with the 0(3) NLSM. This will form the topic of Section 3. 

The Landau-Ginzburg description is an approximate field theory of the magnons and their 
interactions. We represent the magnon field as m. It is a vector under the 0(3) symmetry 
akin to original field n and shares all of n's original discrete symmetries. However m is not 
constrained to live on the unit sphere. The magnons have a gap and a relativistic dispersion 
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relation. The simplest effective Hamiltonian with these characteristics is 

1 

n = -(n2 + (a^.m)2) + — |mp + A|m|^-i7- (mx n), (7) 

where 11 is the momenta conjugate to m. We have added a |m|^ term to ensure overall stability. 
The corresponding Lagrangian (with H in the z direction) is 

1 

C = -((9tm)^ - (9^m)^) + if(ma;(9tmj^ - mj^^tm^) + —(m^ + m^) - —|m|^ - A|m|^. (8) 

The last three terms form the effective potential for the model. 

When H < A, the minimum of the Landau- Ginzburg potential occurs for m = 0. But when 
H > A the minima of the potential now occur for field configurations of the form 

2 2 H^-A^ 

= 0. (9) 

As we are interested in the low-energy behavior of the theory, we focus upon the low-energy 
fluctuations about these field configurations. Introducing the new fields m and $ via 

(m + m,)e±^* = ^^^^^|^, (10) 

with ml = {H^ - A2)/(8A), the effective Lagrangian can be rewritten (to quadratic order in the 
fields) 

C = I((9,m,)2-(a.m,)2) + ((aim)2 + m2(9i<l>)2)-((a.m)2 + m2(9.$)2) 

tt2 

+ 2H{ml + 2mom)dt^ - 2{H'^ - A^)m'^ - -^^l- (^1) 

The low-energy physics in C is governed solely by the field $. The m and modes are massive 
and may be integrated out. Doing so results in a Lagrangian of the form 
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^^[iB,^Y-^mYy (12) 



the bosonic form of the Luttinger model. Here K and vp are the effective Luttinger parameter 
and Fermi velocity. These parameters can be related to the various parameters appearing in the 
original Landau- Ginzburg Hamiltonian. However this can be done only at the mean-field level; 
quantum fluctuations strongly renormalize their values. The approach we thus take is to access 
K and vp directly from the integrability of the 0(3) NLSM. This, as has been indicated, is done 
in Section 3. 

To complete this section we identify the fields of the 0(3) NLSM in terms of $. We can 
factorize the boson, $, into right and left moving pieces each describing the excitations at the 
two Fermi points, $ = 0/? + 0l- The (Euclidean) propagators of 0l/_r are given in terms of K. 
We have 

(0l(2;)0l(w)) = -^\0g{z-w\ 
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{4>r{z)^r{w)) = -^\og{z -w), (13) 

where z/z = vpT ± ix. The fields, n±{x, r), creating 5*^ = 1 excitations are given in terms of the 
boson, $, by 

n± = e±'*, (14) 
and so are governed by the propagator of the form 



{n+{z)n_{w)) = {x + VpT^) 4k. (15) 

We thus understand the spectral functions of low-energy excitations in terms of K when H > A. 
Moreover we can express the z-component of the magnetization operators in terms of the U{1) 
boson currents (in Euclidean space). Fluctuations in the magnetization density are given by 

= = -^9,$ = -id.^n + dM = 7^Ur + Jl), (16) 

while the corresponding spin current, js, equals 

M',^^ = --5.$ = - — {d,c^R - dM = ^{Jr-Jl). (17) 

Current- current correlators then can be deduced from 

{jl{z)jl{w)) = - ^ (18) 
[z — wY 

an immediate consequence of (fTSj) . 

The above identification of the fields is minimal in nature. It is based upon both the Landau- 
Ginsburg analysis together with the simplest association of the fields in the NLSM with those of 
the Luttinger liquid (fT2|) . However it misses subleading terms. To identify such terms directly 
from the Landau-Ginsburg analysis would seem difficult. To circumvent this difficulty, we analyze 
a spin-1/2 ladder system. With ferromagnetic interactions across an individual rung of the ladder 
we obtain an effective spin 1. And with antiferromagnetic interactions along the length of the 
ladder, we obtain an effective antiferromagnetic Heisenberg spin-1 chain. In a magnetic field, the 
ladder system can be reduced to an equivalent spin-1/2 chain. Bosonizing such a chain in the 
standard fashion yields a Luttinger liquid, precisely as in ()12j). Proceeding in this fashion offers 
us the advantage that it allows a more complete identification of the fields, including subleading 
terms. We so find 

K 

M,{x,t) = dr^(x,T) + ccos(2tt(Q(x,t) - Mx)); 

n-{x,T) = e-**(^'")(l + Cicos(27r(e(a;,r) -Mx)). (19) 
Here B(x, r) is the boson dual to that of $ and may be written as 

e(x,r) = -(0i(^)-0R(z)), 

TT 
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while c and ci are some numerical constants. We see that the leading terms agree with that 
of the Landau- Ginsburg analysis. The subleading terms are interesting in that they depend 
upon the magnetization, M, of the system, and so lead to incommensurate behaviour. This 
incommensurate behaviour has been argued to occur [^21 in the context of spin-1/2 ladders with 
ferromagnetic rungs. 

The map from spin-1/2 ladders to effective spin-1/2 chains has been studied extensively. 
Two general approaches have been employed. In the case where the spin-1/2 ladder is strongly 
antiferromagnetically coupled across the rung (we are however interested in ferromagnetic cou- 
plings), a straightforward map takes the ladder in a magnetic field to a spin-1/2 chain [HOI I31j. 
In the case where the spin-1/2 ladder is weakly coupled across the rungs (either ferromagnet- 
ically or antiferromagnetically) more detailed studies have been carried out, using successive 
bosonizations/fermionizations jHSl- Due to the simplicity of the former, we adapt this method 
to ferromagnetic rung couplings, as discussed in Appendix C. 

The need to adapt the map is predicated upon two pieces of physics. The Luttinger parameter 
of the magnetically induced gapless phase of spin-1/2 ladders with ferromagnetic/antiferromagnetic 
rungs satisfies K > 1/ K < 1. However a naive application of the map yields K < 1 regardless 
of the type of rung coupling. Additionally, the map as it stands does not take into account the 
different types of gapless behaviour. The perpendicular susceptibility of ladders with antiferro- 
magnetic rungs, for example, sees power law like decay at /c = (tt, tt) and exponential decay at 
k = (vr,0). But the situation is reversed for ladders with ferromagnetic rungs. 

3 Computing K and from the Integrability of the 0(3) 
NLSM 

In this section we present exact results for the H > He phase of the 0(3) NLSM. They confirm 
the Luttinger behavior implied by the Landau- Ginzburg description discussed in the last section. 
In particular we show how to extract the phenomeno logical parameters, K and vp, appearing in 
this description from the integrability of the 0(3) NLSM. 

We begin by describing how the ground state of the 0(3) NLSM is altered in the presence of 
a magnetic field. Neglecting interactions, the energy of an excitation with bare energy, in a 
finite magnetic field, if, along the z axis is 

E - S,H. 

This leads to a splitting of the spin-1 triplet of magnons. In a mean-field 'semi-conductor' 
picture, this splitting is sketched in Figure 1. As the magnetic field, is increased beyond the 
critical value of H = = A, the ground state of the spin chain dramatically changes. The 
ground state now begins to fill in with a finite density of magnons (and so is akin to a doped 
semiconductor). If the Sz = 1 magnons were non- interacting, they would condense in the lowest 
possible energy level. But the magnons interact repulsively, and so they fill the ground state as 
if they were fermions or hard-core bosons. The technical manifestation of this behavior is the 
minus sign in the S'-matrix at zero momentum transfer (^ = in (j2(J|) ). With H larger than A, 
the energy of an excitation is potentially negative and excitations carrying 3^ = 1 appear in the 
ground state. We illustrate this in Figure 1 schematically by plotting the dispersion relations 
for the three types of excitations in the system: Sz = —1, 0, 1. 
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Figure 1: A sketch of the ground state of the 0(3) non-hnear sigma model in the presence of a 
magnetic field that exceeds the gap, A. The Sz = 1 band is now partially filled with excitations; 
the 5*^ = band is unchanged as it does not couple to the magnetic field; and the Sz = —1 band 
is shifted upwards in energy. 



The ground state of the system for H > He is therefore a sea of magnons. The low-energy spin 
excitations consist of slight deformations of the filled sea in the Sz = 1 band. These low-lying 
excitations can be of arbitrarily small energy and so are gapless. If we linearize such excitations 
above the Fermi energy of this band, we obtain a Luttinger liquid characterized by a Luttinger 
parameter, K, and a Fermi velocity, vp. 

To describe the ground state exactly we need to take into account interactions of the particles. 
These interactions are parameterized by their scattering. As the 0(3) sigma model is integrable, 
the scattering matrix is known exactly. Exploiting the theory's relativistic invariance, we can 
parameterize a particle's energy and momentum in terms of its rapidity 6, defined as 

E = Acosh(^); P = Asinh(^). 

In this parameterization, Lorentz boosts send 6 — » 6 + a. Lorentz invariant quantities like S 
matrix elements therefore depend only on differences of rapidities. The S matrix for scattering 
the two-particle state \ai{9i)a2{92)) (01,02 = 1,2,3) into the final state 1 04 (6'2) 03(^1)) is ^ 

^a3a4^l(^12) + f^aiaa (^0204^2(^12) ~\~ ^aia4,^a2as, cTsi^u), (20) 

where 6 = 61 — 62 and 

27^16 



{6 + m){e - i2Tiy 
2m{i7i - 6) 

[9 + iTi){e - i2Tiy 



Since the ground state is filled with solely 5*2 = 1 particles, we need the S matrix element, 5*++, 
for scattering of two such particles: 

S^M = (21) 
9 + tir 
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This is found from the change of basis, A± = ■^{Ai + iA2), where here Ai is an operator creating 
an excitation carrying quantum number i. 

Knowing 5* allows us to determine the density of states per unit length, p{9) , of the sea of 
particles at if > He. The derivation is standard and can be found for the case at hand in [HH] . 
For completeness we repeat these arguments in Appendix A. At zero temperature, the repulsive 
nature of the particles leads them to fill the sea up to some Fermi momentum, so that the density 
satisfies p{6) = for \6\ > Op- For \9\ < Op, we have 

p{e) = A cosh(^) + r d9'p{9')T^+{9 - 9'), (22) 

where 

T^40)^^^de\ogS^49) = -^. (23) 



The first term on the right-hand-side of (j22j) is the free term, while the kernel, r++, appearing in 
the second term measures the strength of the interactions. As r_|__|_ > 0, the interactions lead to 
a density of states greater than the bare value of ^ cosh(6') associated with purely free fermions. 
The strength of the magnetic field, H, determines the Fermi rapidity, 9f- To find it, it is first 
convenient to introduce the dressed energy, e(^), of the magnons. This is the amount of energy 
the system loses when a particle of rapidity, 9, is removed from the sea. It is given by 

e{9) = Acosh{9)-H+ / d9'ei9')T++{9 - 9'). (24) 

At zero temperature, the particles fill all allowed levels np to 9 = 9^. Therefore the Fermi 
rapidity is determined by solving (j^^ subject to the boundary condition, 

e{9F) = 0. 

Thus for 16*1 < 9f, ^{9) < 0. The energy of the system per unit length at zero temperature is 
given by 

E{H) = / p{9) [A cosh(0) -H]=A — e{9) cosh(0). (25) 



J -Op J -Bp 27r 

The zero-temperature magnetization follows immediately from the above equation, 

M{H) = -dnEiH), (26) 

with the corresponding susceptibility given by x{H) = OhM^H). 

We now are able to compute the Luttinger parameters K and vp appearing in the effective 
theory as functions of H/A. The key is to study this system in terms of its gapless quasi-particle 
excitations present when H > A. These excitations are not the original magnons, but rather the 
excitations above the magnon sea. Since they are gapless, they move either to the right or left at 
the Fermi velocity. We define their rapidity 9 via e(^) = ±vfp{9) = /le^^, where the mass scale 
p is arbitrary and can be redefined by a shift in 9. To compute vp for the gapless excitations 
above the magnon sea, we first note that 

_ de _ 09 de 
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Figure 2: A plot of the Fermi velocity of the low-energy excitations present in the Sz = 1 band 
when H > A. 

Here p is the dressed momentum of the excitations and is given by an equation similar to the 
one governing the energy: 



where (j){6) = logS^^{6)/2'Ki is the unrenormalized scattering phase. By comparing (f^Hj) with 
p2|l we see dgp = 27rp. The Fermi velocity is then given by 



While the integral equations and cannot be solved in closed form, it is easy to solve 
numerically. This yields the curve plotted in Figure 2. Its most notable feature is the vanishing 
of vp for H near A. Below we will give a power-series expansion for vp, valid when H is near A. 

The coupling K is similarly straightforward to find. As we now show, it is related to the 
renormalized spin, Spi, of excitations near the Fermi surface. Sr measures how interactions 
change the response of a spin 3^ = 1 excitation near the Fermi surface to a change in the 
magnetic field. It is related to the dressed energy by 




(28) 



Vf = 



1 de 
2np{e) de 



e=eF 



(29) 



Sr 



(30) 



and so obeys the integral equation. 




(31) 
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Figure 3: A plot of the Luttinger parameter describing the low energy excitations present in the 
Sz = I band when H > A. 

The H dependence of Sr arises from the dependence of 6f on H. 

The basic idea behind the relationship of K and Sr is that the low-energy excitations near 
9 = Op are free-fermionic. That is, their dressed S-matrices are —1 [S^. In contrast to 
here it is relatively easy to demonstrate this and we do so in Appendix B. As a consequence 
of the free-fermionic behavior, the only non-vanishing matrix elements of the current operator, 
ji = —2K(vp^dr — idr^)(j)L , involve a single particle-hole pair. Lorentz invariance requires that 
these matrix elements be given by 



where 6^. and Op are the rapidities of the particle and hole respectively. The constant c is simply 
related to both K and Sr. In a Luttinger liquid, K appears in the current- current correlator 
(fTHj) . To relate this to c, we insert a complete set of states and use the matrix element (P^. 
giving 



On the other hand, c is related to Sr, as integrating the above matrix element gives the value 
of Sr on the one-particle states: 



{OhOp\jLm^) = l^ce 



{dn+ep)/2 



(32) 




(33) 
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—I 



dx 



27r 

-zc6{ei-e2). 



^^g(e2+9i)/2giMx(e«2-e9i) 



Thus c = i27iSii and so 



K = Sl 



(34) 

Given 5*7? is described by the integral equation (|^. i^T follows immediately. The results are 
plotted in Figure 3. 

It is straightforward to find explicit power-series expansions for K and vp valid for H near A. 
We begin by computing e{6). For H slightly larger than A, we can expand e{6) around 6 = 0: 



eie) = do + d2e^ + oie^). 

Plugging this into (j21I) demands do and d2 satisfy 



(35) 



do 
d2 



{A-H){1 + 2T++{0)9f) + A 



A 



(36) 



where r_|_+(0) = l/vr^. To determine 6p we apply the condition e{6p) = 0, resulting in 



er = a/2(| - 1) + ^r+4o)(f - 1) + o((f - if') 



Similarly we can show p{d) to be 



PO + P20'^ 



(37) 



(3^ 



where 



27rpo = A(l + 2r++(0)^f + 4^^r^+(O)) + C»(0:^); 

We can find the renormalized spin similarly. 

These expansions give the power-series expansions for K and fj? to be 



K 



Vf 



A 



V2 



1/2 



3tt^ \A 
8 ^H_^ 



■fj \ 3/2 



37r2 V A 



3/2 



(39) 



(40) 



We have restored the bare spin- wave velocity vfq (earlier set to one). 



13 



We see that as a function of the magnetic field, the Luttinger parameter is 1 at threshold 
H = A, and then increases with increasing field strength. That ir=latif = Aisa universal 
result already obtained by j^HlllZI- It is independent of the exact nature of the excitations that 
fill the ground state at finite magnetization. In the low magnetization regime, the excitations 
of the ground state interact only weakly due to their low density. Thus to a zeroth order 
approximation, they can be treated as free fermions with a corresponding Luttinger parameter 
of K = 1. Thus in the analysis of antiferromagnetic spin ladders -ft" = 1 at if = A is again 
found even though K decreases (opposite to the behaviour of the spin-1 chains) with subsequent 
increases in H. Similarly in mapping the spin chain onto an interacting Bose gas [34j, one arrives 
at = 1 for if ~ A. Here however K at least increases with increasing field. By fine tuning 
the strength of the interactions in the gas, one could in principle come close to reproducing the 
behaviour as predicted by the 0(3) NLSM. 

We illustrate this in more detail. The equations governing the energy and renormalized 
spin of the Bose gas are |35| (in the case of the Bose gas the quantity equivalent to Sr is the 
renormalized charge) 



Here c is the strength of the interactions in the Bose gas. If c = tt, we see the equations of the Bose 
gas are nearly identical to those of the 0(3) NLSM. The sole difference (apart from some trivial 
shifts and rescalings) lies in the dependence of the bare energy upon 6, i.e. eQ°^'^{6) = 6^ — H 
while e^^^^(^) = Acosh(0) — H. At low energies the two are identical, but diverge at higher 
scales. 

4 Scaling Behavior at > A 

Using the results for K and vp derived in the last section, we describe the scaling behavior 
present in the magnetized phase of the spin chain. 

4.1 Susceptibility and Spin Conductance 

We first consider the computation of the magnetic susceptibility. To take into account the effects 
of a magnetic field we add a term to the Lagrangian of the form 



r++(^^) 




(41) 



TT C2 + ^2 • 




(42) 
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Here we have used that the oscillating term of vanishes with the integration over x. Doing 
so the relevant correlator in computing the susceptibility is then 

This correlator can be evaluated in the bosonic formulation of the Luttinger liquid giving 



ttvf Av 



1 1 



(M,(x,r)--a,$(0)) = — ^ -TTTl^ + -TYTl^ ■ (43) 



sinh^C^z) sinh2(^^)J 



Here we have evaluated it at finite temperature. By taking the appropriate analytical continu- 
ations, we can arrive at the retarded correlator necessary to compute the susceptibility, x{H). 
We so find 

K K(H) 
ttvf 7rvF{ri). 
We point out that the static susceptibility is independent of temperature. 

As we are able to compute the susceptibility of the spin chain directly from the thermody- 
namics of the 0(3) NLSM, we are able to perform a non-trivial check on the correctness of this 
calculation. The susceptibility, ^it T=0 is given by taking two derivatives of the energy, as 
displayed in and (f^ . Expanding this in a power series as explained at the end of the last 
section gives the energy per unit length to be 



A2 

EiH) = -- 
SO that the susceptibility is 



25/2 X X 2 ^5/2- 



(45) 



xiH) ^ -OlE^H) ^ ^ - + ^ + O - 1 V ^ (46) 



We see that we recover the susceptibility as derived from the above Kubo formula. 

We can also compute the spin conductivity. The spin conductivity measures the resulting 
spin current arising from a spatially varying magnetic field and is defined via 

The spin current operator has the general form 

K 

is = —vf — + term dependent upon M. 

71 

However the latter term in js involving M will not couple to 6C in ()42|1 . As such the spin 
conductance is given simply from a Kubo formula, 

Rea,{Lj,H) = \imjlm{jj-^dt<i'))^^^{u,k) 



VfTT 

k^O TT^ k 



lini —-lm{dx<^dt^)^^,{uj, k) 



= vfK{H)6{uj). (47) 
This result is valid at any temperature within the Luttinger liquid framework. 
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Figure 4: Plots of the scaling exponent rj for the transverse spin-spin correlator, (^^..(a;) 5*^.(0)) ~ 
\x\~^. In the top panel rj is plotted versus the applied magnetic field while in the bottom panel 
the rj is plotted versus the magnetization of the ground state. In this latter graph we also plot the 
values (squares) of 1] found in Campos Venuti et al. through a density matrix renormalization 
group analysis. 

4.2 Inelastic Neutron Scattering 

We can make a series of predictions relevant for inelastic neutron scattering experiments with 
H > He- For scattering near wave- vector k = n, the experiments probe the single particle spec- 
tral weight of the field, n°'{x,t). The scattering cross section is given in terms of the correlation 
function, 

a{iu,k)(x /"dcjrffce-'"*+''="^(n'^(x,t)n'^(0,0)). (48) 

Because 77.3 is massive, its spectral weight is exponentially suppressed at low energies and tem- 
peratures and for such cases we need only consider the contribution of the correlator {n^n~): 



a{uj, k) (X dxdt e 



-iujt+ikx , 



n+(a;,t)n-(0,0)). 



(49) 



We first examine the T = behavior of this correlator. 
From the relation (fT^ . the correlator {n^n~) equals 



(n+(x,r)n"(0,0)) 



2+^;^r2)V4i^ 



[X 



+ 2cIcos{2txMx)- 



(X2+4T2)^+V4if 



(50) 



The exponent, r] = 1/2K, governing the leading piece of the transverse spin-spin correlator is 
plotted in the top panel of Figure 4 as a function of magnetic field. We see that its H = A+ 
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value is 7] = 1/2, equal to that of free fermions. For values of H slightly larger than A, r/ = 1/2K 
takes the form 

1 1 / 2^/^ ,H .I/O 8 , H . „ , H , q /r,\ 
" = 2]f = 2 - — ( A - 1' ' + 3^< A - 1' + °( A - ^) ' )■ '^1) 

In the bottom panel of Figure 4 we plot rj versus the magnetization of the ground state. We 
also compare our computation to that of numerical simulations of a lattice integer-spin chain, 
H = J J2i Si-Si+i done by Campos Venuti et al. jH]. The agreement between the two is reasonable 
but not overwhelming. At larger values of the magnetization, M, this is to be expected. For a 
lattice spin chain the magnetization must saturate at some critical value of the applied magnetic 
field. He ~ J. In contrast the magnetization of the 0(3) sigma model has no such upper bound. 
Thus we must work at values of the applied field where the magnetization of the system is 
small or alternatively, if -C J. Given that A ~ AJ in a spin 1 chain, this limit can only be 
ambiguously met at best. At smaller values of the magnetization, the disagreement may arise 
from differences between the lattice theory and its continuum version. Nonetheless, we find 
reasonable agreement. 

At finite temperature we can easily determine the form of {n~^n~). Via a conformal trans- 
formation we have, 

^2 1/4K 1 



(n+(x,rK(0,0)) = (^) 



Analytically continuing and then Fourier transforming thus gives us an expression for the cross 
section at finite T, 

1 / TT X-2+1/2X f3 2 1 / ^ y2+K+l/AK 

X (fK+i/,K{-u^, ^{k + 2nM)) + fj,_,y,i,{(luj, ^{k - 27rM))) , (53) 
with f^{x,y) equal to jH], 

h,{x) = Re[(2^)^i?(l^,l-7)], (54) 

where B is the beta function, B{x,y) = r(x)r(y)/r(x + y). We so obtain a scaling form 
for a{LJ,k) with f{x,y) a universal function. The scaling form two parts: one relevant for 
wavevectors k near and one yielding a contribution for k near the incommensurate wave vector 
27rM (understanding that all wavevectors involving have been shifted by vr). 

We can analyze the above expression for a{uj, k) in the small and large T limits. In the case 
T <^ M and uj,k <^ M, we can safely ignore the contribution of the second term in ()53j) . We 
then find 

a{u},k) ~ —^^^^^Q{-u}~VFk)Q{vFk-u)\^ ^ j ; T-^u,k; 
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O / OtT \ -2+1/2/!' TT 1 1 

~ T-^+y^". (55) 

In the first case, T <^ uj,k, the leading term is temperature independent. For the behavior in 
the second case to be observed, we need A ^ T. Otherwise thermal excitations involving the 
other two bands, (5*2 = 0,-1), would alter the Luttinger liquid behavior of the ground state. 
For wavevectors near k ± 27rM the second term in (j53p dominates and we find instead 

a{u;,k±27TM) ~ ^^r2(l/4ir + K) ~ VFk)Qivpk + u) 

X ; T'^uj,k; 
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1 is: . 1 



We thus see the power law dependencies of (y{uj, k) near k ~ 27rM differ from those near A; ~ 0. 

For scattering near the wave- vector, k = 0, the spectral weight of the magnetization operator 
is probed. It is this operator that creates/destroys excitations near A; = 0. At low energies and 
temperatures the cross section, a, is given by 

a{uj,k)oc y"dte-^'^*+^^^(M,(x,t)M,(0,0)). (57) 

The other two spin components of the magnetization are massive and do not contribute at low 
energies. By the fluctuation-dissipation theorem, we can recast a in terms of the imaginary piece 
of the corresponding retarded correlator: 

^) °^ - ^_g-/,^ MM,M,)^^J^, k). (58) 

In the case that uj,k <ti M, we obtain 

K io 

a{uj, k) oc — —^{5{ijj + vpk) + 5{uj — vpk)). (59) 

For wavevectors near 27rM we find instead 



a(uj,k±2TTM) ~ ^— -r'(l - K) sin(7r(K - i)) ^^'^^''^ > ( - 

1 — e~'^'^ vp V 



X (sgn(w + A;) + sgn(w — A;)); T<tiu!,k 
2c^uj ( 27r \2^^-3sin(7rfs:) K 



a{u;,k±2nM) ~ -^^^{j^) ~ '-^^^B\^,1 - K)^{K/2y, T»u;,k;{60) 
where ipi^x) = log(r(x)). 
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4.3 NMR Relaxation Rates 



As with the inelastic neutron scattering cross section, we can make a series of predictions for 
NMR relaxation rates. The NMR relaxation rate is given in terms of the spin-spin correlation 
function [TH]: 

^ —A^p{k)A^^{-k){n''n^)ik,u;N), (61) 



Ti ^ J 2ti^ 

/3=1,2,3 



oc A„/3AaT,(n^n^)(x = 0,u;^ ~ 0). (62) 



where ujn = InH is the nuclear Lamour frequency with '-y^ the nuclear gyromagnetic ratio and 
the Aaf3 are the hyperfine coupling constants. In the above we have assumed H is aligned in 
the 3(z)-direction. We assume that the hyperfine couplings are independent of the wavevector 
k. Hence 

1 

We now consider the correlator {nn) in more detail. Specifically we will compute the contribu- 
tions to T^^ coming from the different branches of low energy excitations. 

If the hyperfine couplings are such that the transverse fields are dominant we obtain for 1/Ti, 

1 1 i^^o J 

= S 2^ j /^i/2i.(/?^/vr) + ^ j /.v2™(/3../vr), (63) 

where h^j^K is as defined in ()54|). In the low temperature limit the second term is subdominant 
and we have 

In the low density limit (i.e. H ~ A+), we have 1/Ti ~ T^^l"^ . In this limit 1/Ti has a 
dependence upon T identical to that of antiferromagnetic ladders [211 E21- However as H is 
increased, the power law dependencies of the two cases diverge. 

If on the other hand the hyperfine couplings are such that the contribution from (M^M^) is 
important we find. 



;i ~ lim /dte^-*(M,(t,x = 0)M,(0)), 



As -f^ > 1 for any finite H in excess of A, the first of the two terms dominates. 



4.4 Validity of Luttinger Liquid Picture 

The Luttinger liquid picture we have developed is precise only at zero temperature. Nonetheless 
we have argued that this behavior will persist to some degree at finite temperature. We are in 
position to analyze qualitatively at least whether this is indeed true. To this end, we compute the 
susceptibility at finite temperature and finite magnetic field. We do so using a more sophisticated 
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Figure 5: Plots of the susceptibility as a function of applied field for a variety of temperatures. 



formalism than presented previously: the thermodynamic Bethe ansatz. For the 0(3) sigma 
model, the appropriate equations were originally given in [21] but can also be found in [T^ liU].^ 

Plots of the exact susceptibility for various temperatures derived from these equations are 
given in Figure 5. Note in particular how the peak at = A appears as the temperature 
is lowered. At T = 0, we see from (|46|) that this turns into a square-root singularity. The 
divergence of the susceptibility from its T = value is a reasonable indication of the persistence 
of the Luttinger-liquid picture at finite temperature as x can be given directly in terms of K and 
vp (see (j21|))- Thus we expect finite temperatures to be destructive of Luttinger liquids associated 
with fields, if, only slightly in excess of the gap where a finite temperature drastically rounds 
off the square root singularity appearing in x- However at higher fields where the density of 
magnons in the ground state is larger and so presumably more robust against small temperature 
perturbations, the susceptibility is equal to its zero temperature value. In this case we expect 
the T = Luttinger liquid picture to remain valid. 

We note as an aside the square root singularity in x appearing in Figure 5 at T = is a 
generic feature. It appears in any model (not only spin chains) with a gap and a low-energy 

^Apropos of nothing, we have also computed the susceptibihty at large temperatures but zero field. In the 
zero- field large-T limit, Sachdev and Damle provide a high temperature computation of the susceptibility 
predicated on integrating out higher Matsubara frequencies. To subleading order it is given by 

From a numerical analysis of the TBA equations at temperatures in the range, T ^ lO'^A — lO^A, we find 
agreement with this analytical form at the 1%-level. To obtain this agreement it was important to include the 
subleading term. 
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Figure 6: A plot of the dependence of the 3D ordering temperature, Tc/J, as a function of 
apphed field, H, for two different spin chain compounds. The relevant values of J for two 
materials are JcsNiCh ~ l^ii' and JAgVP2S6 ~ 320K. The corresponding values of inter to 
intrachain anisotropics are J' /J ~ .017 (CsNiCl^) and J' /J ~ 10^^ {AgVP2SQ). The left 
ordinate scale corresponds to AgVP2SQ while the right scale corresponds to CsNiCl^. 



0.04 




quadratic dispersion relation. Its appearance is not related in any way to the 0(3) NLSM 
being integrable (although the integrability allows us to compute exactly the coefficient of the 
divergence). We thus expect this square- root divergence to appear also in spin chains with large 
easy-axis anisotropics. Indeed the rounded finite temperature counterpart of this divergence has 
already been seen in NENP |37j . 

Beyond the effects of finite temperature, the Luttinger liquid phase can be destroyed by the 
presence of intrachain couplings. In mean field theory, three dimensional order due to such 
couplings develops at a transition temperature, Tc, given by 

l = |aJ'x±(g = 0,^ = 0,T = T,)|, (67) 

where x± is the staggered susceptibility of a single chain, 

X±= f Me*'^*-*'"(n+(x,t)r2-(0))retarded, (68) 



and J' is the intrachain coupling. Here we explicitly include the lattice cutoff, a ^ ~ J. g is the 
deviation from the wavevector, vr. From the expressions for given in (|19|). we obtain, similar 

to US], 



B{ t , 1 )B{ t , 1 

^8K 47rT ' 4:K' ^8K A-kT ' AK' 



J. I 



27raw,2+2X+^)) Sin(7r(i^+ 47^)) \ K ^ ^ _ „. ^f(^- - 27rM) ^ _ ^ _ ^ 



B(— + . - 

^ 2 Sir AtiT ' AK' 
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-4 -s^-'^^^^'-^-i^) <-) 

where u± = {uj/vp^k). Combining this expression with the mean condition leads to a transition 
temperature given by 

Tc_ ^ 1 Vp ^J' VfQ ^■■^2K/(4K~1). 

S = sin(— — ,1 - — ). (70) 

In deriving this expression we have ignored the term in x± appearing in ()7Up dependent upon 
M whose contribution is neghgible for small T. 

Plotted in Figure 6 is the dependence of Tc upon H for the spin chain compounds CsNiCl^ 
and AgVP2SQ. We see that as H is increased there is a corresponding increase in T^. Thus at a 
fixed temperature and sufficiently large interchain coupling, J', we expect that as H is increased 
there will be a transition to a state with long range order. We thus see that the Luttinger liquid 
phase is stable at a given temperature for smaller fields. 

As discussed in the introduction, field induced long range order has been observed in the 
spin chain material, NDMAP IHS]. Although this compound has a large single ion anisotropy, 
we make an attempt to compare experimental observations with our analysis. From inelastic 
neutron scattering experiments, the spin coupling, J, of NDMAP is J = 2&.bK [Hj. For a 
field applied along the chain (perpendicular to the easy plane defined by the anisotropy), the 
corresponding gap A is given in terms of J and D = .3 J, the strength of the single ion anisotropy, 
by 03], 

A ~ .4J - MD = 5.3K = 3.75T, 

with the Lande g factor equal to g = 2.1. In order to obtain a reasonable fit with the experimental 
observations we take J'/ J = 2 x 10"'^. The results are plotted in Figure 7. The value used for J'/ J 
is considerably smaller (by a factor of 6) than the value employed in These authors, however, 
were able to match the experimental data by assuming x± took the form x± ~ -yyj^j^ 
A some field independent constant. We however find that A is both field dependent and generally 
larger in magnitude than the value of A used in jHS]. To compensate in matching the data, we 
need to take J'/ J to be smaller in magnitude. 

At fields not far in excess of the gap, we find reasonable agreement between the data and 
the theoretical computation. At larger values of the field, H, the computed Tc exceeded the 
observed value. This is not surprising. As we saw in our computation of the exponent t], the 
assumptions underlying our computation fail to accurately describe the physics at large fields. 
With H/A > 2, we enter a region where the magnetization per site of the system is some 
significant fraction of its maximal value of 1, a regime our method, which puts no bound upon 
the magnetization, is unable to handle properly. 

The destruction of the Luttinger liquid phase due to an applied magnetic field occurs not 
only for spin-1 chains but for antiferromagnetic spin-1/2 ladders [SI]. Unlike the result found in 
the case of two antiferromagnetically coupled spin-1/2 chains j3T], we find however is strictly 
monotonic as a function of applied field. The difference in the two cases is again that we have 
K > 1 and increasing with field while in the case of antiferromagnetic spin ladders, K < 1 
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Figure 7: A plot of the dependence of the 3D ordering temperature, Tc/J, as a function of 
apphed field, H, for the spin chain compound, NDMAP. 



and decreasing with field. There is an additional instability that affects antiferromagnetic spin 
ladders, an instability to the formation of a spin-Peierls state due to coupling to the lattice jlB] . 
However spin-1 chains are thought to be robust against this perturbation 46j. 
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A Analysis of the Ground State at T=0 



In this section we determine the characteristics of the zero temperature ground state in a mag- 
netic field, H, exceeding the gap. At the heart of this analysis lies the insistence that the ground 
state wave function be compatible with the scattering amplitude, 5++, as defined in (jTIjl . This 
amplitude characterizes the wave function, iIj{6i, . . . ,6^), of a ground state with N particles 
through the application of periodic boundary conditions. In allowing the i-th particle, 6i, to 
traverse the entire length, L, of the system and so commute with the other — 1 particles, the 
wave function picks up a phase 

N 

The first term is a product of scattering amplitudes arising from the i-th particle scattering with 
the remaining — 1 particles. The second term comes from the bare momentum carried by the 
i-th particle. With periodic boundary conditions, we must have 



Taking logarithms of this constraint leads to the quantization condition, 

LA 



sinh(e,) + 5^-^log5++(%), (72) 



2n V ^ 2m 

3+^ 



where A^j is the quantum number characterizing the rapidity, 6'j, (i.e. log(l) = i27cNi). As 
we approach the continuum limit, we can then write a difference equation between successive 
integers, Ni and A'j+i = A'i -|- 1, 

1 = ^cosh{d,)6d, + V -L^delogS++){0,,)6ef, 66, = 9,+, - 6,. (73) 

The density of states per unit length at 6, is defined as n{6i) = 1/L66i. With this we can rewrite 
the above in the continuum limit as, 



(74) 



where p is the density of occupied states and p is the density of unoccupied states, so that 
p + p = n. 

At zero temperature, the equation for p can be simplified. With T = 0, the ground state has 
a Fermi surface at 6'i?. Excitations with rapidities, 16*1 > Op, do not appear in the ground state. 
Thus 



n{e), \e\<eF; 

0, 1^1 > 9f; 
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[n{e}, \e\ > Of; 

and we obtain equation for p at zero temperature. 

We have computed the interacting density of states. We can also compute the interacting 
energy of the 5*^ = 1 excitations in the ground state. The total energy of the system equals 

E = j d9{Acosh{e) - H)p{9). (76) 

If we vary the density of occupied particles and holes, p ^ p + 6p and p ^ p + 6p, the total 
energy varies accordingly 

5E= [ de{Acosh{9) - H)5p{0). 



But we can also express this variation in energy in terms of the interacting or dressed energies 
of the excitations: 

5E = j de{5p{e)t+{d)-5p{B)t-{e)). 

e^{9)/t^{9) mark the energies needed to excite a particle/hole above the ground state. At zero 
temperature they are defined such that 



>0, 


1^1 


> Of 


= 0, 


1^1 


< Of 


= 0, 


1^1 


> Of 


<0, 


1^1 


< Of 



Together make up a monotonic, smooth function via e = + e . Comparing the two 
expressions for SE and using the constraint on the variation in p and p coming from ()74|1 . 



5p + Sp = J d9'5p{9')T++{9-9') 



we arrive at equation 



B Determination of Dressed Scattering Phase 

In this appendix we compute the dressed two body scattering phase of a particle-hole excitation. 
To compute this phase we examine how the momentum of one of the added excitations is altered 
by the addition of both excitations to the system. 

The two-body scattering phase of the particle-hole excitation is defined by 

Sphi9p, 9n) = L{ppi9p) - ppo(^p)), (78) 

where Pp{9) is the momentum of the particle when the particle/hole pair is present, while Ppo{9) 
is what the momentum of the particle would be without the hole and without the effect of the 
particle's presence on the sea of excitations already present in the T = ground state. 
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Let 6p, 6h be the rapidities of the particle and hole we intend to add to the system. Without 
having added either excitation, the dressed momentum of the particle is given by 

PpoiOp) = Asinh(^p) + 27r / d9no{9)4>{9p ~ 9); 

m = ^ log S+^{9), (79) 



where uq = p + p is the density of states given in (j22|) . This is the density of states without 
either the particle or hole present. Once the particle and hole are added, the density of the 
ground-state sea is altered. Rather than being given by (j221), it is governed by 



n{9) = ^ cosh(5) + ^(r++(^ - 9p) - S{9 - ^,)) + / d9'n{9')T^+{9 - 9'). (80) 



'F 



The 1/L term (r(^ — 9p) — 5{9 — 9h)) represents the disturbance the particle-hole excitation 
produces in the density of states. The 5{9 — 9h) arises from removing the particle to create the 
hole, while T{9 — 9p) arises from including the particle in the integral j d9 pTj^j^: the presence 
of this particle alters the density of states in the same way the continuum of particles in the 
ground state does. 

Taking into account the alteration the addition the particle-hole excitation has upon the 
density of states, the momentum of the added particle becomes 



Pp{9p) = Asmh{9p) +2n d9n{9)(f){9p - 9). (81) 

J -Op 



"Bp 
-dp 

As is unshifted by the addition of the particle-hole pair we are able to write n{9) as 



n{9) = no{9) + ^n,{9) - ^5{9 - 9^). (82) 

Substitution of (jH^ into (jHn|) then allows the scattering phase to be reduced to 

fdp 

Sphi9p,9h) = -2n(l){9p-9h) + 2Tr d9n,{9)(l){9p - 9); 

J -dp 

n,i9) = T^+{9-9p)-T++{9-9j,)+ r d9'n,{9')T++{9-9'). (83) 

J-ep 

We are interested in computing the scattering phase right at the Fermi surface, i.e. 5ph{9p = 
9f, 9h = 9p). But we then immediately see 

V(^f,^f) = -2710(0) = -7r. (84) 

Hence the low-energy particle-hole S-matrix is, 5* = e*'^'''"^^^'^^-' = —1, as we claimed. Here 
we have focused upon the changes in the momentum of the particle. If we had examined the 
momentum of the hole instead we would have arrived at an identical conclusion. 
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Figure 8: A two rung spin ladder. 

C Reduction of a Ferromagnetic Spin- 1/2 Ladder to an 
Effective Spin-1/2 Chain 

In this section we consider a map reducing a spin-1/2 ladder with ferromagnetic rung interactions 
in a magnetic field to an effective spin-1/2 chain. The essential idea behind the map is already 
discussed in a number of articles (SUl l^- However there the treatment considered a ladder 
with antiferromagnetic rung interactions. Here we adapt the map to ladders with ferromagnetic 
rungs. 

For the sake of simplicity we consider a two-rung ladder, as pictured in Figure 8. The 
Hamiltonian of the ladder is given by, 

Ti = Hj + Hjj^ 

= J{Sal ■ Sa2 + Shi ■ Sh2) 

+H{S:i + Sli + S:2 + Sl2). (85) 

We will suppose that J > while J^±_ < J±z < 0. 

We first treat the limit J = where the rungs decouple. The eigenstates on each rung with 
their corresponding energies are then 

1+) = ITT); 



4 



H 



|0) = _(|ti) + |iT)); 



J 



±± 



ii); 



J 



±z 



H 



\s) = -^(ITi)-liT)); 



J 



±_L 



4 ' 



^6) 



where the first three states are that of the triplet while the remaining state correspond to the 
singlet. The two states |+) and |0) are of lowest energy. It is these two states that we will use 
in constructing the effective spin-1/2 chain. 
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To take into account Tij we perform second order perturbation theory. We find that the 
energies of the four basis states describing the low energy Hilbert space of both rungs, | + +), 
1/72(1 + 0) ± |0+)), and |00), are shifted as follows: 



I + +): 
1/72(1 + 0) ± |0+)) : ± 



J 
2 

J 
2 

^ 2(J_L_L — Jj^z) 8J_L_L ^ ^ 

Understanding |+) = ||) and |— ) = ||), we can write down an effective system of interaction 
spin- 1/2 's: 

^e// = Y^^t^2 + 4^r) + J>l^l - H\al + a'^) + c, (88) 
where the various effective couplings are given in terms of the original couplings by 

Jl = J; 

Jl = J/2 + a; 

We point out that as a can be negative, so can J^. negative is essential to correctly describing 
the physics of the spin-1 chain. 

Having reduced the spin ladder with ferromagnetic rung couplings to a spin-1/2 chain we 
now bosonize. Applying the standard results (see [2]) we find that the associated Lagrangian is 



K 

2^' 



K = (90) 

while the fields are given by 

n-{x) = ^(5-(a;) + 5fe-(x)) =a-(x) =e-**(^)(l + cicos(27r(e(x)-Ma;))); 
v2 

M,{x) = S'^{x) + S'^{x) =a'{x) + l/2 = — —dr^x)+ccos{2n{e{x)-Mx)). (91) 

VpTl 

Here M is the magnetization density of the spin-1/2 ladder (and so the corresponding spin-1 
chain). In terms of the magnetization, M^, of the effective spin-1/2 chain, M = + 1/2. We 
see n~ is defined in terms of the ky = component of the ladder spin operators. Here we thus 
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expect to encode power law correlations for wavevectors near k = (7r,0). In a ladder with 
antiferromagnetic rungs, power law correlations would be seen instead for k = (vr,7r). 

With < (provided a is sufficiently large) we see that K > 1 in correspondence with both 
our analysis and the results of |33i- This marks out the main difference between ladders with 
ferromagnetic rungs and ladders with antiferromagnetic rungs. With antiferromagnetic rungs, a 
similar analysis leads to > and so < 1 (|30, (STj). 
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